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1.  Summary 


The  first  part  of  the  research  is  that  we  have  expanded  the  Exact  Scientific  Computational  Library  (ESCL), 
and  the  Dixon’s  algorithm  [I]  on  rational  N  by  N  matrix  inverse  was  implemented.  We  studied  the  relation 
of  required  length  A/ofp-adic  expansion  and  the  prime  p  by  experiments,  and  the  possible  use  of  the  length 
of  periodicity  of  a  rational  number’s  p-adic  expansion  in  determining  the  length  of  required  M  in  rational 
matrix  operations. 

The  second  part  of  the  work  is  to  develop  and  implement  computational  algorithms  for  p-adic  cyclic  code 
generation,  which  is  based  on  the  results  of  the  paper,  Modular  andp-adic  cyclic  codes,  by  A.R.Calderbank 
and  N.J.A.  Sloane  [2],  Algorithms  and  software  have  been  developed  to  give  an  alternative  solution  to 
factorize  the  polynomial  X"-l  over  the  ring  of  integers  modulo  pa,  where  p  is  a  prime  not  dividing  n,  and  it 
can  generate  the  table  of  cyclic  codes  using  the  divisors  of  )C-I  as  their  generator  polynomials. 

All  the  implementation  of  ESCL  is  in  C++,  as  well  as  the  software  to  generate  p-adic  cyclic  codes. 

2.  Progress  on  ESCL  Using  P- adic  Arithmetic 

In  the  design  of  ESCL  we  employ  several  developing  tools  and  a  specific  library--NTL.  We  use  Visual  C++ 
6.0  and  MFC  to  develop  and  debug  the  codes.  MATLAB  Symbolic  is  used  for  comparison.  We  adopt  the 
NTL  library  to  process  the  input  and  output  of  arbitrary  long  integers,  which  are  the  numerator  and 
denominator  of  a  rational  NTL,  written  by  Shoup  [3],  a  high-performance,  portable  C++  library  for  number 
theory  that  provides  both  data  structures  and  algorithms  for  arbitrary  length  integers.  NTL  allows 
manipulation  of  integers  for  vectors,  matrices,  and  polynomials  over  finite  fields,  and  arbitrary  precision 
floating-point  arithmetic* 

Our  ESCL  has  the  following  operations; 

1)  Converting  rational  to  p-adic  number  and  via  verse. 

2)  Single  rational  number’s  calculation:  addition,  subtraction,  multiplication  and  division. 

3)  The  addition,  subtraction  and  multiplication  of  two  rational  matrixes, 

4)  The  inverse  of  a  single  rational  N  by  N  matrix. 

In  order  to  improve  the  efficiency  of  ESCL,  a  proper  prime  p  and  the  length  M  ofp-adic  expansion  need  to 
be  determined.  We  did  experiments  on  a  random  rational  matrix  with  different  prime  numbers  and  record 
the  execution  time.  Our  experiments  show  that  the  running  time  is  more  related  to  the  length  Mofp- adic 
sequence,  where  “AT  can  be  calculated  based  on  the  prime  p: 

m  =  2[log(&T' )/  log  p\,  (l) 

where  the  S  =  f" [  X,  and  Xt  is  the  Euclidean  length  of  the  ith  column. 

The  process  of  choosing/?  and  M  can  be  done  as  follows.  Suppose  that  for  certain  "p'\  we  can  get  “AT  very 
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smalt.  If  we  choose  a  small  Mt  then  use  ( 1)  to  calculate  prime  p.  Thus  we  will  get  a  relatively  big  prime  and 
a  theoretic  smallest  M,  and  some  adjustment  is  needed  during  the  process.  The  flow  chart  is  shown  in 
Figure  l. 


Figure  I-  Estimating  the  prime  number 

Dixon's  algorithm  for  matrix  inverse  from  uExact  Solution  of  Linear  Equations  Using  P-adic  Expansions" 
[1]  was  implemented,  which  is  the  newly  added  function  since  the  last  report.  To  get  the  A'\  for  which 
AjcA  J~/f  we  take  the  following  three  steps, 

1)  Get  the  nxn  matrix  Ct  whose  entries  lie  in  [0,  p-I] 

AA‘i=I*>  AA'1  mod  p  =  I  mod  p=>AC  =  I  mod  p;  =>  C  -  A'1  mod  p. 

2)  Compute  a  p-adic  approximat  ion  x1  to  A'1,  x1  is  a  nxl  column  of  A'1  . 

flt-1  JW-I 

x'=YsXiP'  =^Cb‘  m0d  P^P' 

j-0  1-0 

3)  Transfer  x*  to  rational  field. 


Repeat  doing  2)  3),  until  bi=bn-]t 


T 

'o' 

"O' 

and  b0= 

0 

b,= 

1 

...  bn.i= 

0 

_o_ 

_0_ 

_1 

3 


Figure  2  shows  an  example  for  this  operation. 


Please  enter  tlie  operation  numher<such  as  1  or  2>:  10 
Please  input  the  division  of  the  s square  matrix: 


Please  input  ROD  i  elements : 
L,5,2 


Please  input  ROU  2  elements: 

Li, 7 

Please  input  ROU  3  elements: 
3,8,2 

[he  prime  number  is:  5 
Result : 

[9/8, -1/8,-11/16 
L/24, -1/24, 5/48 
-1/6,1/6,1/12 1 


Figure  2,  Matrix  inverse  example 

Our  system  is  capable  of  carrying  out  exact  arithmetic  operations  for  rational  numbers  using  the  finite 
segmented  p-adic  numbers.  The  representation  of  the  rational  numbers  has  several  appealing  properties. 
Given  the  usual  representation  of  the  positive  integers,  it  is,  in  a  genuine  sense,  the  natural  extension  to  the 
rational.  The  algorithms  for  addition,  subtraction,  and  multiplication  are  those  of  the  usual  integer 
arithmetic;  the  division  algorithm  is  truly  the  analogue  of  multiplication. 

When  two  p-adic  rational  numbers  are  added,  subtracted,  multiplied,  or  divided,  the  result  is  an  infinite  but 
eventually  repeating  sequence  of  digits  (periodicity). 

Despite  the  advantages  of  finite  segment  p-adic  arithmetic  and  exact  representation,  there  are  certain  issues 
need  to  be  addressed: 

1)  Improve  the  EstimatcMQ.  Currently,  the  A/ estimation  does  not  always  work  well.  In  some  situation,  it 
cannot  give  the  sufficient  number  of  digits  for  p-a die  sequence  for  exact  computation.  To  avoid  this 
kind  of  errors,  we  have  increased  its  size  to  30  and  more. 

2)  A  more  efficient  method  for  prime  number  selection  for  the  given  rational. 

3.  The  Length  of  Periodicity  for  the  P-adic  Expansion 

To  determine  what  is  the  sufficient  length  AY  of  p-adic  expansion  for  a  rational  number  in  the  matrix 
operation  is  not  a  trivial  problem.  Let  us  observe  what  happens  after  the  arithmetic  operations  of  two  p-adic 
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sequences. 


All  rational  numbers  can  be  uniquely  written  in  the  form  of  a  =  . 

We  know  that  a  real  number  is  rational  if  and  only  if  its  decimal  expansion  is  periodic.  Similarly,  a  p- adic 
number  is  rational  if  and  only  if  its  /7-adic  expansion  is  periodic.  Consequently,  since  we  are  primarily 
interested  in  the  p-adic  expansions  of  rational  numbers  we  will  be  dealing  only  with  /vadic  expansions 
which  are  periodic. 

The  expansion  eventually  repeats  to  the  right.  That  is,  if  a  is  a  rational  number,  then  it  has  a  repeating 
pattern  of  djPJ  in  its  P-adic  expansion,  Le.,  it  is  of  the  form 

a=,Aft...Asa0...a^i 
For  example: 
i  =  .231(p=5) 

|  =  .4l3(p=5) 

The  operation  of  addition,  subtraction,  multiplication  and  division  in  the  set  of  p-adic  numbers  are  quite 
similar  to  the  corresponding  operations  in  decimals.  The  main  difference,  however,  is  that  we  proceed  from 
“left  to  right"  rather  than  from  “right  to  left”  as  we  do  with  decimals. 

Addition/subtraction 

Assume  that  we  have  two  P-adic  sequences:  (s<t) 


a  =  ,At 

...Asa,...a 

b  =  A 

Line  up  these  two  sequences 

AXAV. 

■Ard\ .  ..cjj-j 

4 

■A 

Set  a,_,+l=Ct 

‘  •  ^l  +  ft  =  Cl 
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cr..cncv..c„... 


A^...Axav.n,_t 
B,  ... B, 


VAW 


Let’s  consider  the  right  side  of  the  vertical  stroke  line 

C[  C2  C3  ...  Cj|  c,  .  cn... 

+  b|  b2  b3  ...  . .  bm  b|  ...  bm... 


C|+b|  c2+b2  c3+b3  ...  Ck+bi,... 

Suppose  two  integer  x,y  satisfy  the  condition  that  x  ^  y  and  cx+bx  =  cy+by 
V  c*  and  by  can  be  any  number 

in  the  worst  case,  we  must  make  sure  that  c ,  =  ay  =  a ,  and  cx  =  by  =  b, 
x  =  kJ*n  +  i=  k2xm  +j 
y  =  Uxn  +  /=  k4  x  m  +  j 
=>  x-y  =  (kl‘k3)  x  n  =  (k2-k4)  x  m 
*.■  x  *  y 

kl-k3  and  k2-k4  are  nonzero  integers 

x-y  must  be  exactly  divisible  by  n  and  m,  the  smallest  integer  which  satisfies  this  requirement  is  LCM(n.m), 

The  carry  digits1  modification 

We  only  need  to  consider  the  (t+l)th  cany  digit,  because  it  will  cause  the  exception  of  the  periodic  part. 
Since  the  maximum  carry  digit  of  the  addition  of  two  numbers  is  1,  we  can  classify  the  following  numbers 
into  three  types  and  they  have  different  effects. 

Type  One;  The  (t+l)th  digit  is  less  than  p-l,  assume  the  (t+l)th  digit  is  q: 

\'q<p-\-=>q  +  l<p 

=>  .q...  + 1  (the  carry  digit)  =  .{q  + 1)... 
so  only  one  digit  will  be  affected. 

Type  Two:  The  next  r  digits  following  the  t th  digit  have  the  value  p-l 
Assume  the  (t+r+ 1  )lh  digit  is  q,  and  q  <  p  —  1 

_ r _  r 

{p  - !)...(/?  - 1  )q...  + 1 . . .  =  .0...0(g  + 1). . . 

r  + 1  digits  will  be  affected. 

In  addition,  r  <  LCM(n,m) 


Because  if  r=  LCM(n,  m) ,  then  it  belongs  to  the  type  three. 

Type  Three:  The  digits  following  the  \(h  digit  are  rotation  numbers  of  p-l 


6 


.p  - 1  + 1  =  .0  =  0 

The  result  is  0. 


We  can  draw  a  conclusion  that  the  maximum  length  of  the  p-adic  expansion  is: 

2  x  LCM  («,  m)  +  max(5 ,  t)  - 1 . 

Multiplication 


axb  =  .Al...Asal...anx.Bx...Blbl...bm 


=  .At...AIai:..anx.Bl...B,+.Ai... A, a,  ...a„x .0...06, ...bK 


=  .Al...AIai...a„x.Bl...B,  +.Ai...Asx.0...0b]...bm+.0...0ai...a„x,Q...0bi...bm 

^ _ _ _ v - - - v - “ — *  v - — v - # 


Let's  analyze  the  three  parts  separately, 
Part  1: 


Z 


***  — As  .  *  *  dn  x  *  *  -  Bt  —  .Cj  * . .  Cu  Cj . , .  cn 

That  is,  the  result  is  a /?-adic  sequence,  which  has  a  periodic  part  of  n  digits. 


Part  2: 

v  ...Bf  e  Z 


.Af . .. As  x .0„.06t ...b„=.Df..Dvdx...dm 

That  is,  the  result  is  a  p-adic  sequence,  which  has  a  periodic  part  of  m  digits. 


Part  1  +  Part  2: 


.At . ..  Asa, ..  ,a„  x  .5, . ..  B,  +.At . . .  A,  x  .0...06,  ...bm 

' - * - y  ' - v - - ' 

1  2 

According  to  the  previous  conclusion  we  get  from  the  Addition  part. 


=>  x.B,  „aBt  k .0, .-0^  *  +  -  .E{...Ewe j 

This  is  a  j^adie  sequence  which  has  a  periodic  part  of  LCMfn.m)  digits. 


Part  3: 
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...a„  x.0.,.06,  ...bm 

=  (.a,.,.an)x(l  +  p”  +  p2n  +  ...  + p'°)x(jbi...bm)x(\  + pm  +  p2m  +...  +  p°  )x/rs+' 
=  (.al...a„)x(.bl...bm)x  p’*‘  x(l  +  pn  +  p2n  +  ...  +  p“)x(l  +  pm  +  p2m  +  ...  +  /?") 

t ^ - *  * - — - V- - - - 

4  5 


Part  4  {.a]..,an)x(J>r..bm)x  p'"  eZ . 


Part  5  is  the  one  we  need  to  focus  on. 


0  +  A 

'  +  p2n+, 

..  +  p*‘)x(\  +  pm  +  p 

2m  + 

.  +  #>•) 

n- 

■]  nTt 

m- 

1  n-1 

*5.10.. 

.010.. .0... 

X.10...010...0. 

- 

a\ 

<*2 

<*1 

«4 

«s 

«6 

a ,  a* 

ak... 

xft, 

A 

A 

A 

A 

A 

bj  bK 

A- 

a\h\ 

<*A 

aA 

<*A 

aA 

<*A 

aib\  a%bx 

aA-- 

aAt 

aA 

aAi 

asb2 

Qffii  Qjbj 

ak- 1 A  ■  ■ 

aA 

aA 

a,  by 

aAs 

ak~Ai  •  •  * 

«A- 

Ignore  the  carry,  the  k//i  products  digit  is  akbt  +  oi_[b2  + . . .  +  Q\bk 

Set  (k'+l)x  LCM(n,m)>  k  >  k'xLCM{n,m)  let  LCM(n,m)  =  £ , 

akb\  +  ak_t  b2  +  ...  +  ak_t+lbt  +  ak_tbl+ 1  +  ok_t_A+  2  + . ..  +  ak-it+Au  + ...  +  Qk-kAk-t+x  +  ■  +  aiA 

s — - - v — - '  ^ - v - — — '  r  ' 

\  2  3 

**'  Qk +1  =  ak-t+i  —  =  Qk-k'i+t  ^  bf  =  =  ”■  ”  ^k-f+t 

.■.part  lean  be  rewrite  as  ak_k,b\  +ok_k.l_lb1  +  ...  +  aA-k-t  +  afbk-k  7+,  + ...  +  ak_t.Mb, 
part  2  can  also  be  rewrite  as  ak_t.f bx  +  ak.k-t-Ai  +  ...  +  £*A-*7  +afbk_k.f+ ,  +  ...  +  a1_*.,+,6, 
part  3  can  be  rewrite  as  ak^k.,b\  +  ak_k  f_lb2  + , . .  +  atbk_k.t 
We  can  find  that  part  1  equal  to  part  2,  and  because  there  are  LCM(n.m)  items  in  this  part,  we  must  get 

a,  X  ft,  once,  which  is  the  only  item  has  the  value  1,  and  all  other  items’  value  is  0,  that  is  the  result  of  pan 

1  is  1,  so  far  as  part  2. 

Part  3  is  the  actually  the  front  part  of  part  t ,  that  means  the  value  of  part  3  may  be  1  or  0. 

Base  on  this  discovery,  we  can  divide  the  configuration  \n\o  LCM (n, m)  x  LCM (n,m)  modules.  Let’s 

take  100100100... x  101010, ..(p-  5)  as  an  example. 
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LCM(2,3)=6,  so  we  divide  the  module  in  6  x  6 
.100100100100100100100100100100100100 


r* 

.10101010. 

1OIOI01010101010101010101010 

1  0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

I 

0 

0 

] 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

l 

0 

0 

\ 

0 

0 

1 

0 

0 

I 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

J 

0 

0 

] 

0 

0 

1 

0 

0 

1 

0 

0 

I 

0 

0 

1 

0 

0 

I 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

I 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

t 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 

0 

] 

0 

0 

1 

0 

0 

I 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 
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4 

4 
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5 

4 

5 

5 

5 

5 

We  can  find  two  type  of  modules  here,  one  is: 

10  0  10  0 

0  0  0  0  0 

10  0  1 

0  0  0 

1  0 
0 
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This  one  is  corresponding  to  the  part3; 

The  other  one  is: 

10  0  10  0 

0  0  0  0  0  0 

0  0  1  0  0  1 

0  0  0  0  0  0 

0  10  0  10 

0  0  0  0  0  0 

This  one  is  corresponding  to  the  parti  part2  and  so  on,  we  can  accumulate  every  line  of  this  module,  the 
result  is  always  I ,  which  is  the  previous  conclusion  we  draw  for  part  l . 

Ignore  the  carry,  assume  that  ck  is  the  kth  digit  of  the  product,  we  can  find  that  ci+,  —  ck  +1 


Now  let’s  continue  to  get  the  result: 

101111  212222  323333  434444  545555  656666  767777  878888  989999... 

101 1 1 1  212222  323333  434444  001111  212222  323333  434444  001111... 

100100100... xioioio... 

=10111 12122223233334344440 

Conclusion:  In  multiplication,  the  length  of  periodic  part  of  the  product  is: 

LCM(m,n)x(p  - 1) ,  (2) 

where  m  and  n  are  the  length  of  periodic  part  of  the  two  multipliers,  p  is  the  prime. 

The  length  of  periodicity  of  the  resulting  p-adic  sequence  can  be  very  large  from  (2).  But  if  we  should 
represent  all  the  p-adic  sequences  with  a  complete  period  during  all  the  calculations,  we  will  definitely 
carry  out  all  the  arithmetic  operations  exactly. 

4.  Computational  Algorithm  and  Software  for  Generating  P-adic  Cyclic  Codes 

A.R.Calderbank  and  N.J.A.  Sloane  [2]  "Modular  and  p-adic  cyclic  codes  ",  studied  how  to  lift  from  binary 
cyclic  Hamming  code  of  length  7  to  octacode,  which  is  equivalent  to  the  binary  nonlinear 

Mordstrom-Robinson  code,  and  provided  a  general  structure  of  cyclic  codes  over  Z6,  Z\b . then  to  the 

2-adic  integers  Z:°o.  That  is  to  say,  all  the  cyclic  codes  share  a  common  generator  polynomial  X1  +  x:  + 
(Jt-l)x  -1,  which  satisfies  X1  -  X  +  2  =  0,  X  is  a  2-adic  number. 

The  purpose  of  our  study  is  to  provide  an  alternative  solution  to  generate  the  p-adic  cyclic  codes,  which  can 
be  implemented  efficiently  using  computer  science  methodology,  and  we  developed  software  to  generate 
p-adic  cyclic  codes. 

Theoretical  Background 
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According  to  the  definition  of  cyclic  code,  a  subset  S  of  Zqn  is  cyclic  if  (an.|,ao,ai . £^2)6  S  whenever 

(a0,a,!...,afl.2,an.|)  e  S.  A  linear  code  C  is  called  a  cyclic  code  if  C  is  a  cyclic  set. 

We  also  know  that,  R  =  Zq[X]/(Xn- 1 )  forms  a  polynomial  ring.  By  defining  the  following  map  between  Zqn 
and  Z,[X]/(X"-!), 

it:  Zq"  Z„[X]/(XM),  (ao.a, . a,.,)  |->  ao+a^x+..,+afl.1xn•,  (3) 

We  can  set  up  a  one  to  one  correspondence  between  an  element  of  space  Zqn  and  an  element  polynomial  of 
ring  Zq[X]/{XM). 

Connect  ideals  of  ring  Zq[X]/(Xn-l)  and  cyclic  codes  contained  in  Zqn  by  the  following  theorem; 

Theorem  1  Let  ti  be  the  linear  map  defined  in  (3).  Then  a  nonempty  subset  C  of  Zqn  is  a  cyclic  code  if  and 
only  if  n(C)  is  an  ideal  of  Zq[X1/(XM).  (Proof  in  [4],  p. 1 36). 

For  the  convenience  of  our  discussion,  we  give  the  following  definition  of  the  generator  polynomial  of  a 
cyclic  code. 

Definition  1  The  unique  monic  polynomial  of  the  least  degree  of  a  nonzero  ideal  /  of  Zq[X]/(Xn-l )  is  called 
the  generator  polynomial  of  /.  For  a  cyclic  code  C,  the  generator  polynomial  of  tt(C)  is  also  called  the 
generator  polynomial  of  C 

By  using  this  definition,  we  can  set  up  an  one-to-one  correspondence  between  the  cyclic  codes  in  Zq"  and 
the  monic  divisors  of  XM  e  Zq[x], 

To  be  precise,  each  monic  divisor  of  is  a  generator  polynomial  of  a  cyclic  code  in  Zq 

For  example,  by  factorizing  the  polynomial  x6  -  l  e  Z2[x]: 

X6-  1  =(!  +x)2(l  +  x  +  x2)2 

We  can  list  all  the  monic  divisors  of  x6  -  1 : 

1,  I+x,  l+x  +  x\ 

(1+x)2,  (I  +x)(l  +  x  +  x2),  (1  +  x)2(l  +X  +  X2), 

(1+x  +  x2)2,  (1 +x)(l +X  +  X2)2,  l+x6. 

These  nine  monic  divisors  are  related  to  nine  ideals  of  ring  Z}[x]/^-l,  and  thus  related  to  nine  cyclic  codes 
of  length  6.  Based  on  the  map  ?i,  we  can  generate  all  these  cyclic  codes. 

For  instance,  we  can  get  the  cyclic  code  corresponding  to  the  polynomial  (I  +  x  +  x')'  as  its  generator  by 
the  following  way: 

Step  I.  Expand  the  generator  polynomial,  (1  +  x  +  x2)2=  I  +  x2  +  x4. 

Step  2.  Multiply  the  generator  polynomial  with  all  the  polynomials  within  the  ring  Z^xJ/X6-!,  and  get  the 
following  ideal  /. 

{0,  I  +  x2  +  X4,  X  +  x3  +  Xs,  I  +  X  +  x2  +  x3  +  X4  +  Xs}. 

Step  3.  Using  map  n,  we  get  the  corresponding  cyclic  code: 

{000000,  101010,010101,  1 1 1 1 11}. 
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In  general,  we  denote  [n.k,dj  to  describe  a  linear  code  and  naturally  a  cyclic  code,  where  n  is  the  length  of 
the  codeword,  k  is  the  number  of  the  base  vectors,  and  d  is  the  least  distance  between  any  two  code  words. 

The  following  theorem  can  relate  [n.k.d]  of  a  cyclic  code  to  the  features  of  the  generator  polynomial  of  an 
ideal  I  of  ring  Zq[X]/(X-i). 

Theorem  2  Letgfx)  be  the  generator  polynomial  of  an  ideal  of  2q[X])(X'-l).  Then  the  corresponding  cyclic 
code  has  dimension  k  if  the  degree  of  g(x)  is  n-k,  and  the  length  of  the  codeword  is  n. 

In  the  above  example,  the  generator  polynomial  is  g(x)  =  /  +  +  x4,  and  n  =  6,  k  =  6-4  =2. 

This  relates  to  a  [6.2,3]  cyclic  code.  The  length  of  the  codeword  is  6,  and  this  cyclic  code  will  have  2k  =  22 
=  4  code  words. 

The  following  theorem  decides  the  number  of  cyclic  codes: 

Theorem  3  Let  x"  - 1  e  Zq[x]  have  the  factorization 

r 

xB- 1  =nPjd<x), 

i=l 

where  pi(x).  p2(x), ....  p/x)  are  distinct  monic  irreducible  polynomials  and  e,  £  1  for  all  i  =  1.2, . r.  Then 

there  are  17,  ,]  (e,+i)  cyclic  codes  of  length  n  overZ,. 

The  key  point  of  finding  all  the  cyclic  codes  over  Z”  is  the  factorization  of  X'-i  over  Zq.  We  will  discuss 
our  method  to  factorize  the  polynomial  X'-l  over  ZH. 

Mechanism  to  Generate />-adic  Cyclic  Codes 

The  idea  of  generating  p- adic  cyclic  codes  is  given  by  the  following  theorem  : 

Theorem  4  Let  q  =  pa.  I  £a<oo.  lfg,(x)  e  Zp[x]  is  a  monic  irreducible  divisor  of x”-l  over  Zp,  then  there 
is  a  unique  monic  irreducible  polynomial  g„(x )  e  Zq[x]  which  divides  xn-l  over  Zq  and  is  congruent  to  gi(x) 
mod  p. 

The  complete  proof  can  be  found  in  [2],  If  we  define  a  cyclic  code  of  length  n  over  Zr ,  we  can  get  another 
cyclic  code  of  Zpl  whose  generator  polynomial  is  obtained  by  lifting  the  generator  polynomial  of  the  first 
cyclic  code  to  Zp2 .  And  then  we  can  continue  this  way  to  get  the  cyclic  code  over  where  !  <a  <x 

An  Alternative  Algorithm  for  the  Generation  ofp-adic  Cyclic  Codes 
Our  algorithm  has  two  steps: 

a)  Factorization  of  X'-l  over  Z2; 

b)  Lifting  the  generator  polynomial  over  Z2r  and  then  to  Z2r+i. 

Step  a)  is  the  foundation  of  our  implementation.  In  order  to  simplify  our  situation,  we  will  let  n  be  a  prime, 
which  has  the  form  8m- 1,  so  that  X-l  can  always  be  factorized  over  Z2  in  the  following  form; 

xn-l  =(x-  l)xt(x)ft2(x) 

where  all  the  factors  are  irreducible. 
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Factorization  of  A"-/  over  Z2 
Since 

xn-l  =(x-  l)ni(x)n2(x) 
jii(x)  =  ao  +  a]X  +  a2x2  +  ...  +  an.^x"'1'2 
Jt2(x)  =  b0  +  b|X  +  b2x2  +  ...  +  bn.t/ix"‘l,: 


a„  b,  e  Z2  ( 0  <  i  <  n-I/2  ), 

So 

xn-1  =(x-  1)(  I  +  x  +  x2  +  ...  +  x"'’)  =  =(x-  1)JI,(x)7I2(x) 


Let 


aob0  =  1  (mod  2), 
aob|  +  a|b0  =  I  (mod  2), 


aob„.i/2 +  aib(n-i/i>.)  +  ...  +  a^i^bo  -  1  (mod  2), 
a|b„.i/2  +  ajb(n.|/2).i  +  ...  +  an-i/2b i  =  1  (mod  2), 


an.i/jbn.,/2  =  1  (mod  2). 

And  the  vector  (aO,al,..„an-l/2,bO,bl,...,bn-l/2)  which  satisfies  all  the  equations  above  will  be  the 
coefficients  of  irreducible  polynomials  Jtt(x)  and  x2(x). 

Therefore,  it  is  easy  for  computers  to  enumerate  from  (0,0,0,...,0)  to  (1,1,1 . 1);  and  because  x"-l  has  the 

unique  factorization,  so  the  computer  can  stop  wherever  it  finds  a  solution. 

Lifting  the  Generator  Polynomial  over  Z2r  to  Z^r+l 

From  theorem  1,  we  know  that  for  every  monic  irreducible  divisor  h,(x)  of  A"-/  over  Zp,  there  is  a  unique 
monic  irreducible  polynomial  in  Z^xJ  which  divides  A"-/  over  Zp,  and  is  congruent  to  h,(x)  mod  p.  Let 

x"-l  =  (x-  l)7i,(x)7t3(x); 

Assume 

ttj(x)  “  f2 1  (x),  7t2(x)  =  g2!(x); 


Let 


f:2(x)  =  f2l(x)  +  2*{c0  +  ctx  +  ...  +  cn.IQxn'l/2), 


gj2(x)  =  g2l(x)  +  2*(do  +  d,x  +  ...  +  d„.l/2xn'ia). 


0  <  ci,  di  <  1, 0  5  i  £  n*l/2. 


We  have 


f22(x)  =  f2 1  (x)  (mod  2),  g>2(x)  a  g2l(x)  (mod  2). 


Let 
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fj2(x)  =  c’d  +  c’,x  +  ...  +  c,n.iaXB'l,21 
gj2(x)  =  d’0  +  d']X  +  ...  +  d’n.i/jx'"1 2. 

Again,  we  can  use  computer  to  enumerate  vector  (Co,C|,...,cn.|/2,do,d[,...,dn.|,i)  from  (0,0,...,0)  to  {1.1,...,  I  >  to 
find  a  value  which  can  satisfy  the  following  equation  group: 

c’nd’o  =  1  (mod  22), 
c’od’i  +  c’|d’o=  I  (mod  22), 

t-  od’n-t/2 +  c'id’(„.|/2).|  +  ...  +  c’n.i.id’g  =  1  (mod  2"), 
c’id’n.i.2  +  C':d'(n.|/2>|  +  ...  +  cVt/jd’i  =  1  (mod  2~), 

cVi/2d5n-u2=  1  (mod  2!). 

The  solution  will  be  the  unique  factorization  of  x"-I  over  2‘,  and  after  getting/:. if.*)  and  g}2(x ),  we  can 
continue  this  step  to  get  fja(x)  and  g2a(x),  (  2  <  a  <  o°  ). 

Software 

The  following  are  some  screenshots,  which  can  demonstrate  our  programs: 


Figure  4  is  the  screenshot  when  p  =  2.  n  =  7,  r  =  8;  p  should  be  a  prime  number,  n  should  be  a  prime 
number  which  has  the  form  8m- 1,  r  could  be  any  integer  larger  than  zero.  We  can  find  that  it  has  the  same 
result  as  in  A.R.Calderbank's  paper. 
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c*  C:\Docuftien1s  and  SetMng3\himter\My  Document&Whiifll  StiiUin  ProJectsVlhesisVleteaseMhe.,,  HE3B 


Command >>1 

Input  prime  number  p  ae>  2 

Input  the  dei/rett  of  polynomial  XAn-l  V 

Input  tl m  powef  to  uhich  the  polynomial  XAn-l  he  lifted  ”>  8 

jp  -  2P  n  *  7,  i4  3  8 

Got  one  solution!!? 

The  value  of  xl  M  I  mod  2  ist  1  +  1ka2  +  lxA3 
The  value  of  *2141  mod  2  is:  1  +  ixAl  +  lxA3 


Cot  one  solution!?! 

[I he  value  of  x!141  mod  4  is; 

[The  value  of  x2141  nod  4  is: 

Col  one  solution?!? 

The  value  of  xll41  nod  8  is: 

The  value  of  x2!4)  mod  8  is: 


-  l  *  2xa1  -  lxA2  +  lxA3 
*  1  #  lxAl  *  2xA2  *  lxA3 


-  1  *  2xAl  +  3xA2  +  lxA3 

-  1  -  3xAi  -  2xA2  +  ixA3 


k]ot  one  solution! tt 

The  value  of  xlMl  nod  18  is:  -  1  —  8x  l  -*  5x  2  +  lx  3 
The  value  of  x2141  nod  16  is :  -  1  *  5xAl  *  6xA2  +  lxA3 


Got  one  solution!!? 

The  value  of  xil4l  nod  32  is: 

The  value  of  x2f4l  nod  32  is: 

Got  one  solution!?! 

The  value  of  xl  14 1  mod  64  is: 

The  value  of  x2!41  mod  64  is: 

Got  one  solution!!? 

The  value  of  xl(41  nod  128  is 

Tlie  value  of  x2  ( 4 1  mod  128  is 

Got  one  solution!!! 

The  value  of  xl 141  nod  256  is 

[f he  value  of  x2l41  nod  256  is 


1  -  6xAl  -  SxA2  +  ixA3 
1  *  5xAl  *  6xA2  +  ixA3 


1  +  26>T1  +  27xA2  +  lxA3 
t  -  27xAl  -  26xA2  +  lxA3 


1  -  38xAl  -  37*A2  +  lxA3 
1  +  37xAl  +  38xA2  +  lxA3 


1  *  9llxAl  *  91xA2  *  lxA3 
1  -  Vlx^i  -  90x~2  ♦  lxA3 


Figure  4.  Output  When  p  =  2,  n  =  7,  r  =  8 
Figure  5  is  the  screenshot  when  p  =  3,  n  =  l  J,  r  =  8: 


Figure  6  is  the  screenshot  when  p  =  2,  n  =  23,  r  -  8: 
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a  C:\Documents  and  SettiitRsVJUiritef  Wy  OocomentsWKual  Studfo  Pro  jet:  is  \1hesis\Re  lease  Vi  be%  H.exe 


The  main  function  of  our  software  is  the  factorization  of  polynomial  ofX'-l.  From  the  coding  theory,  we 
know  that  cyclic  codes  have  many  good  characteristics  for  error-correction  and  cryptography.  It  may  have 
some  merit  to  investigate  the  p- adic  cyclic  codes.  In  this  study,  by  setting  up  a  one-to-one  correspondence 
between  cyclic  codes  and  ideals  of  a  polynomials  ring,  we  can  use  algebraic  methods  to  explore  the 
features  of  /?-adic  cyclic  codes. 

In  the  future,  we  may  use  this  software  as  a  tool  to  study  the  features  of p- adic  cyclic  codes,  and/or  use  the 
software  to  auto-generate  cyclic  codes  in  any  application. 
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